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Effective  Wake  :  Theory  and  Experiment 
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ABSTRACT 

An  improved  theoretical  method  is  presented  for  com¬ 
puting  the  effective  wake  of  propulsors  operating  in  thick 
stern  boundary  layers  on  axisymmetric  bodies.  The  hydro- 
dynamic  intern'  ti  n  between  the  nominal  velocities  upstream 
of  the  propulsor  and  at  the  propulsor  location  is  assumed  to 
be  inviscid  in  nature  and  the  total  energy  is  assumed  to  be 
conserved  along  a  given  streamline  with  and  without  the  pro¬ 
pulsor  in  operation.  Theoretical  predictions  using  the  method 
are  compared  with  experimental  data  obtained  in  the  United 
States  and  Japan  for  five  different  propulsor/axisymmetric 
body  configurations.  For  all  five  cases  examined,  the  com 
puted  total  velocity  profiles  immediately  upstream  of  the  pro¬ 
pulsor  (with  the  propulsor  in  operation)  arc  in  good  agree¬ 
ment  with  the  measured  values.  In  addition,  the  volume-mean 
values  of  effective  velocity  profiles  computed  from  the 
measured  nominal  velocity  profiles  are  in  good  agreement 
with  the  measured  values  of  the  Taylor  wake  fraction  (l-wT) 
for  ali  five  nominal  wake  distributions  over  a  wide  range  of 
propul sor  thrust  loading  coefficients. 

1.  INTRODUCTION 

The  velocity  profile  at  the  location  of  the  propulsion 
device  in  the  absence  of  a  propulsor  is  called  the  nominal 
velocity  profile.  The  effective  inflow  velocity  distribution  ex¬ 
perienced  by  the  propulsor  depends  on  the  mutual  interaction 
of  the  propulsor  (loading  distribution  and  geometric  charac¬ 
teristics)  and  the  stem  boundary  layer.  It  can  be  significantly 
different  from  the  nominal  velocity  distribution.  In  the  design 
of  a  wake-adapted  propulsor,  the  radial  distribution  of  effec¬ 
tive  inflow  velocity  is  often  estimated  by  ratioing  the 
measured  nominal  circumferential-mean  axial  profile  up  or 
down  by  a  constant  factor.  The  factor  is  sometimes  taken  to 
be  (l-wTV(l-wN),  where  wT  is  the  Taylor  wake  fraction  and  wN 
is  the  measured  volume-mean  nominal  wake.  Naval  architects 
derive  the  Taylor  wake  fraction  from  propeller  open  water 
curves  on  the  basis  of  thrust  identity  between  propeller 
powering  experiments  in  open  water  and  behind  tne  ship 
model. 

The  constant-factor  empirical  approach  for  obtaining 
effective  inflow  distribution  is  not  based  on  a  rational  hydro¬ 


dynamic  theory.  Until  recently,  detailed  velocity  surveys  in 
the  presence  of  an  operating  propulsor  were  not  available, 
and  the  actual  distribution  of  effective  velocity  into  a  pro¬ 
pulsor  had  not  been  examined  fully.  Providing  the  correct 
distribution  of  effective  inflow  for  wake-adapted  propulsor 
design  is  essential  to  meeting  the  ever  increasing  demand  for 
improving  propulsion  performance  and  energy  conservation. 

In  order  for  a  propulsor  to  produce  a  required  thrust  to 
power  a  ship  with  minimum  power  and  minimum  cavitation 
at  a  prescribed  propulsor  rotational  speed,  the  effective 
velocity  distribution  used  in  the  propulsor  design  must  be 
very  accurate. 

In  1976,  a  Laser  Doppler  Velocimeter  (LDV)  was  suc¬ 
cessfully  used  by  Huang  et  al.1  to  measure  velocity  profiles 
very  close  to  the  propulsor.  The  measured  velocity  profiles, 
stern  pressure  distributions,  and  stern  shear  stress  distribu¬ 
tions  with  and  without  an  operating  propulsor  provided  the 
necessary  clues  to  the  proper  understanding  of  the  interac¬ 
tion  between  a  propulsor  and  stern  boundary  layer  on  axisym¬ 
metric  bodies.  The  influence  of  propulsors  on  axisymmetric 
stern  boundary  layers  was  found  to  be  contained  within  a 
limited  region  extending  two  propulsor  diameters  upstream  of 
the  propulsor.  The  propulsor/stern  boundary  layer  interaction 
was  found  to  be  inviscid  in  nature.  The  inviscid  approxima 
tion  computer  program  developed  by  Huang  et  al.1  predicted 
very  well  the  measured  total  velocities  in  front  of  the  oper¬ 
ating  propulsor.  Subsequent  detailed  measurements  of  the 
velocity  profiles  with  and  without  a  propulsor  operating  in 
two  axisymmetric  wakes  were  made  by  Nagamatsu  and 
Tokunaga.2  An  inviscid  approximation  again  predicted  well 
the  measured  total  velocities  in  front  of  the  operating  pro¬ 
pulsor. 

Schetz  and  Favin3  4  have  formulated  a  numerical  pro¬ 
cedure  based  on  the  full  Navier-Stokes  equations  to  compute 
the  flow  near  body/propulsor  configurations.  The  computed 
axial  velocities  at  two  propulsor  diameters  downstream  of  the 
propulsor  compared  satisfactorily  with  the  measured  results. 
This  numerical  procedure  has  not  been  applied  to  the  com¬ 
putation  of  effective  wake  in  design  of  wake-adapted  pro¬ 
pulsors. 

The  influence  of  a  stern-mounted  propulsor  on  the 
flow  field  past  bodies  of  revolution  and  a  flat  plate  was 
measured  by  Hucho,5-fi-7  although  no  attempt  was  made  to  ac- 


tually  calculate  the  effective  wake  distribution.  Other  inves¬ 
tigations  related  to  this  subject  were  made  by  Wertbrecht,8 
Hickling,9  Tsakonas  and  Jacobs,10  and  Wald.11  Methods  to 
estimate  effective  wake  were  proposed  by  Raestad,12 
Nagamatsu  and  Sasajima,13  and  Titoff  and  Otlesnov.14  The 
only  known  previous  effort  to  theoretically  address  this  prob¬ 
lem  is  due  to  D.M.  Nelson*  who  developed  a  computer  pro¬ 
gram  for  calculating  the  effective  wake  from  the  measured 
nominal  wake  and  static  pressure  distribution  across  the 
boundary  layer. 

In  the  present  paper,  an  improved  method  for  com¬ 
puting  effective  wake  distribution  from  the  measured  nominal 
wake  distribution  is  derived.  Serious  effort  has  been  made  in 
this  work  to  compare  the  theoretical  velocity  predictions  with 
velocity  distributions  measured  by  an  LDV  or  a  five-hole  pitot 
probe  immediately  upstream  of  an  operating  propulsor.  For 
all  five  cases  examined,  the  computet!  total  velocities  up¬ 
stream  of  the  operating  propulsors  are  in  good  agreement 
with  the  measured  values. 

The  computed  volume-mean  effective  wake  distribu¬ 
tions  are  shown  to  compare  favorably  with  the  measured 
Taylor  wake  fractions  derived  from  self-propulsion  experi¬ 
ments  of  five  different  nominal  wakes.  A  considerable  amount 
of  experimental  data  and  relevant  computational  results  are 
tabulated  to  jiermit  the  independent  assessment  of  the  pres¬ 
ent  method  by  other  investigators. 

2.  THEORY 

-.1  l'ropulson Stern  Boundary  Layer  Interaction 


where  ip  is  the  stream  function  for  an  incompressible  axisym- 
metric  flow  defined  by 

I  dtp 


I  dtp 


Since  the  flow  velocities  are  increased  due  to  the  ac¬ 
tion  of  the  propulsor,  stream  surfaces  are  shifted  closer  to 
the  body  surface.  As  shown  in  Figure  1 ,  a  typical  stream  sur¬ 
face  moves  inward  from  r  to  r(,  while  the  resultant  velocity  is 
higher  than  the  nominal  velocity.  The  resultant  velocity  up,  as 
would  be  measured  in  front  of  the  propulsor,  will  be  called 
the  total  velocity. 


(  Region  of  Propulsor  Influence 
Extends  to  About  2  Dp  i 

(  Upstream  of  the  Propulsor  I 


The  experimental  data  given  in  References  1  and  2 
allnw'one  to  conclude  that  the  influence  of  propulsors  on 
upstream  stern  boundary  layers  is  detectable  only  within  two 
propulsor  diameters  upstream  of  the  propulsor.  Upstream  of 
the  propulsor.  the  mean  circumferential  velocity.  v#,  is  iden¬ 
tically  equal  to  zero  on  an  axisymmetric  body  both  with  and 
without  the  propulsor  in  ofieration.  The  following  assumptions 
are  made  to  derive  a  theoretical  approximation  of  the  hydro- 
dynamic  interaction  between  a  propulsor  and  a  thick  stern 
Imundai y  layer  upstream  of  the  propulsor:  (a)  the  flow  is  ax¬ 
isymmetric  and  the  fluid  is  incompressible;  (h)  the  interaction 
of  propulsor  and  nominal  velocity  profile  is  considered  to  be 
mviscid  in  nature:  thus,  propulsor-induced  viscous  losses  and 
turbulent  Reynolds  stresses  are  neglected:  (o)  the  conven¬ 
tional  boundary-layer  assumption,  3v,V3x  «  3u,/3r  is 
assumed  to  lie  valid  for  the  nominal  boundary  layer  in  the 
absence  of  a  propulsor;  and  (d)  upstream  of  the  propulsor,  no 
energy  is  added  to  the  fluid  by  the  propulsor,  and  the 
propulsor-induced  velocity  field  upstream  of  the  propulsor  is 
irrotatiorial.  The  theoretical  formulation  of  propulsor/stern 
boundary  layer  interaction  subject  to  the  above  assumptions 
has  already  been  given  by  Huang  et  a).,1  and  will  only  be 
briefly  outlined  here. 

The  vector  equation  of  steady  motion  for  an  inviscid 
fluid  is  gr.err  bv  (see,  for  example.  Reference  1.1) 

.  ^  1 

V  v  u>  =  —grad  H,  (1) 

e 

where  V  is  the  fluid  velocity,  w  =  V  x  V  is  the  vorticity  vec¬ 
tor.  and  H  =  p  ♦  l/2p(^  V  )  is  the  total  head,  with  q  mass 
density  and  p  pressure.  For  cylindrical  polar  coordinates  (r,  0, 
x)  with  V  =  (vr,  v»  =  o,  u„),  the  radial  comfionent  of  Equa¬ 
tion  (1)  may  be  written1  as 

/3u,  3vr\  1  3H  dip  ru,  3H 

'  V  3r  3x  /  g  dip  dr  e  dip  ^ 


Figure  1  -  Definition  Sketch  for  Propulsor  Stern 
Boundary  Layer  Interaction 

Since  no  energy  is  added  to  the  upstream  fluid  by  the 
propulsor  and  since  we  assume  no  change  of  viscous  losses 
due  to  propulsor  induction  effects,  the  total  pressure  head 
within  a  given  stream  annulus  upstream  of  the  propulsor  re¬ 
mains  constant  with  the  propulsor  (denoted  by  a  subscript  p) 
and  without  the  propulsor  in  operation;  thus,  upstream  of  the 
propulsor 


1  dH  _  I  (aH\ 

e  dip  e  'dipfp 


or  from  Equation  (2) 

^  /  dux(r)  __  d  vr(r)  \ 
r  r  \  3r  3x  I 

{  du,,(r,,)  ^  3(vr(rp)  +  vpr(rp)) 
r(,  l  3rp  3x 

where  w8  is  the  0-component  of  the  vorticity  vector.  The 
value  of  Wj/r  must  therefore  be  conserved  along  the  stream 
annulus  with  and  without  an  operating  propulsor.  The  radial 
velocity  with  the  propulsor  in  operation  is  assumed  equal  to  vr 
+  vJir,  where  vr  is  the  radial  velocity  without  propulsor  and 
vpr  is  the  circumferential-average  propulsor-induced  radial 
velocity.  Since  the  propulsor  induced  velocity  field  is  assumed 
to  lie  irrotational,  we  require 


<Vx  Vpl)„ 


3x 


3u„ 

drp 


=  0  , 


where  Vp,  is  the  resultant  propulsor-induced  velocity,  and  u, 
is  the  circumferential-average  propulsor-induced  axial  veloci¬ 
ty.  Substitution  into  Equation  (3)  yields  a  simple  relationship 
for  the  location  of  the  new  stream  surface  rp 


•private  rommunioation.  unpublished 


1  3u,  1  /  3  Up  3u»\ 

r  3r  rp  \8rp  3rp  / 

i  .  I  !^1  (4) 

rp  dx  r  dx  1  • 

The  normal  boundary -layer  app-oximation  dv^dx  « 
3u,/3r  is  now  assumed  to  be  valid  for  the  nominal  velocity 
profiles.  In  order  to  check  this  assumption,  the  radial  nominal 
velocity  profiles  were  measured  by  a  TSI,  Inc.  Model  1241 
"X"  -  wire  at  two  axial  planes,  one  immediately  upstream 
(x/Dp  -  -  0.227)  of  the  propulsor  and  one  on  the  propulsor 
plane.  The  values  of  dvjdx  determined  from  velocity  meas¬ 
urements  on  Afterbodies  1  and  2  of  Reference  16  were  found 
to  be  less  than  0.05  3u,/3r.  The  radial  velocity  vr  was  found 
to  vary  very  slowly  across  the  boundary  layer.16  Therefore, 
we  may  assume  that  v^r)  a  v,frp).  Furthermore,  the  value  of 
(r-rp)/rp  represents  the  percentage  change  in  radial  location  of 
a  given  streamline  with  and  without  the  propulsor  operating; 
its  value  is  found  to  be  less  than  0.2  for  most  of  the  cases 
considered.  Thus,  the  order  of  magnitude  of  the  last  term  in 
Equation  (4)  can  be  estimated  as 

I  avr<rP>  _  I  a  Vr<r) 

rp  3x  r  3x 
__ 

r  3r 

and,  therefore  will  be  neglected  in  the  following  analysis. 
Equation  (4)  then  reduces  to 

1  3ux  1  /  3m  3ua  \ 

r  3r  rp  V  3rp  3rp  )  -  ^ 

Within  a  given  stream  annulus  the  massflow  is  constant  with  or 
without  a  propulsor  operating.  Thus,  at  a  given  axial  location  x, 
we  may  write, 


3vr(r)  r-r 


3x 


i 

r  3r 


—  <  0.01 


dw  -  ru,dr  =  rpupdrp.  (6) 

Inserting  this  expression  into  Equation  (5)  yields 

Mu,  =  upd(up  -  uj  .  (7) 

Equations  (6)  and  (7)  are  the  governing  equations  for 
the  propulsor  and  stem  boundary-layer  interaction.  The 
nominal  velocity,  ux(r),  and  the  circumferential-average 
propulsor-induced  axial  velocity,  u,(rp),  can  be  used  to  obtain 
the  new  location  of  the  stream  surface  rp  with  total  velocity 
up.  In  Equation  (6),  the  values  of  r  and  rp  are  equal  on  the 
body  surface.  Far  outside  of  the  boundary  layer,  the  flow  is 
uniform  without  any  vorticity.  Then  we  have  du,  «  0  as  r  -» 

°°  in  Equation  (6),  which  implies  d(Up  -  uj  -  0,  and  up  -  u„ 

-  constant  -  V  (ship  speed),  or  Up  -  V  +  u*,  since  u,  is  zero 
and  Up  has  to  be  equal  to  V  as  r  -* 

Similarly ,  for  an  open  propulsor  in  uniform  flow,  dux 

-  0  in  equation  (7),  and  we  have  Up  =  V  +  u,.  The  total 
velocity  up  for  an  open  propulsor  in  uniform  flow  is  the  sum 
of  the  uniform  ship  speed  V  and  the  propulsor-induced  veloci¬ 
ty  u,  without  any  interaction. 

If  the  nominal  velocity  field  is  irrotational,  the  left- 
hand  side  of  Equation  (3)  must  be  zero;  thus,  the  right-hand 
side  of  Equation  (3)  becomes 

3MP)  _  3vr<rp)  _  gyfrp)  _  0 
3rp  dx  dx 


Using  the  irrotationality  condition  applies  to  the  nominal  veloci¬ 
ty  as  well  as  propulsor-induced  velocity  fields,  e.g. 

aVrp)  _  3u»(rp)  and  3v^rp)  _  3u,(rp) 

8x  3rp  dx  "  3rp  . 

The  above  equation  becomes 

3 

—  [Up(rp)  -  ux(rp)  -  u,(rp)J  -  0  , 
or 

Up(rp)  =  ux(r p)  +  u,(rp)  . 

Thus,  the  total  velocity  up  for  a  propulsor  operating  in  an  ir¬ 
rotational  nominal  velocity  field  is  the  sum  of  the  correspond¬ 
ing  nominal  velocity  and  the  propulsor-induced  velocity  with¬ 
out  any  additional  interaction. 

2.2  Effective  Velocities 


In  current  propulsor  design  and  performance  predic¬ 
tion  practice,  only  the  measured  model  nominal  velocity  pro¬ 
file  is  available.  The  full-scale  nominal  velocity  profile  can 
then  be  estimated  as  the  sum  the  measured  model  nominal 
velocity  profile  plus  an  appropriate  Reynolds-number  correc¬ 
tion.  The  nominal  velocity  profile  for  either  the  model  or  the 
full-scale  propulsor  is  assumed  to  be  available  at  the  outset. 
The  present  theory  can  be  applied  to  obtain  the  effective 
velocity  distribution  using  the  nominal  velocity  profile  as 
input. 

If  the  effective  velocity  is  defined  to  be  the  total 
velocity  with  the  propulsor  in  operation  minus  the  propulsor- 
induced  axial  velocity,  u„  =  u(1  -  ua,  then  Equation  (7) 
becomes 


Mux  =  (u,.  +  ujdu^  .  (K) 

At  large  radial  distances  from  the  propulsor  axis  (r  —  °°),  the 
value  of  the  propulsor-induced  velocity  becomes  zero  and. 
hence,  the  effective  velocity  is  equal  to  the  nominal  velocity  as 
r  —  00 .  The  finite-difference  form  of  Equation  (8)  can  be  written 
as 


|u*. . ,  +  uv,)|u,,.  ,  -  uX|]  = 

K-l  +  Ml  +  Uv,  +  Uajlu,..,  -  u,.j 


19) 


u“. . 


Ua, 


In  practice,  the  nominal  velocities  are  known  from  the  inner¬ 
most  radius  of  the  propulsor  to  a  radial  position  (e.g.,  r  5  3.0 
Rp)  large  enough  for  the  propulsor-induced  axial  velocity  ua  to 
be  zero.  The  value  of  the  effective  velocity  at  that  point  is 
equal  to  the  nominal  velocity.  This  condition  can  le  used  to 
solve  the  effective  velocity  in  Equation  (9)  step  by  step  from 
this  initial  radial  position  inward  towards  the  body  surface. 
The  nominal  velocity  at  the  hub  of  the  propulsor  is  usually 
taken  to  be  the  linearly  extrapolated  value  of  the  measured 
nominal  velocities  near  the  body  surface  rather  than  zero  in 
order  to  avoid  an  unnecessarily  large  number  of  computa¬ 
tional  grid  points  near  the  body  surface. 

Since  the  induced  velocity  us  is  not  known  initially,  the 
present  theory  can  be  applied  in  an  iterative  procedure  to  pro¬ 
pulsor  design  of  performance  prediction.  The  circumferential 
average  propulsor-induced  axial  velocity  u,  will  be  assumed  to 
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be  the  same  within  the  same  stream  annulus  with  and  without 
the  propulsor  in  operation  in  every  iteration. 

In  conventional  propulsor  design,  when  the  nominal 
wake  distribution  is  given  and  either  the  propulsor  thrust  at  a 
given  speed  or  the  delivered  shaft  power  is  given,  the  follow¬ 
ing  procedure  may  be  followed: 

(a)  First,  estimate  the  thrust  loading  coefficient  Cjs  of 
the  propulsor  and  then  estimate  the  propulsor-induced  axial 
velocity  by  using  the  actuator  disk  approximation,17  i.e., 

r 

—  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

R„ 


— —  0.147  0.206  0.254  0.290  0.311  0.315  0.293  0.233  0 
VC'ts 

An  initial  estimate  of  effective  velocities  can  then  be  obtained 
by  solving  Equation  (9). 

(b)  The  conventional  propulsor  lifting-line  computional 
method18 1!,~"  niay  be  used  to  design  a  propulsor  to  operate  in 
t he*  estimated  effective  velocity  profile  u  with  a  prescribed 
nondimensional  circulation  distribution  G  which  produces  the 
required  thrust.  An  improved  solution  estimate  of  the  circum¬ 
ferential-average  propulsor-induced  axial  velocity,  ua,  can  be 
computed  from  the  output  of  the  lifting  line  computer  pro¬ 
gram1"  by  a  simple  integration,-1 


tion  will  be  extrapolated  linearly  toward  the  wall,  resulting  in 
a  nonzero  velocity  at  the  wall.  Since  u,  and  u^  can  be  approx¬ 
imated  locally  by  a  linear  function  of  r  and  rp,  the  mass  flux 
within  the  stream  tube  annulus  given  by  dr  -  r, ,  ]  -  r,  and 
drp  ”  rPi  ♦  l  ~  rt>i  can  be  integrated  from  Equation  (6)  as 


Ii  * 1  In.  i 

Vi.i  ~  Wi  -  AVj  =/  ru„dr  =/  rpUpdrp 

^r,,, 

with  u(r)  =  u,  +  [(r  -  r,)/(r, , ,  -  r,)J  (u,  + ,  -  u,)  and  r,<r<r, . , 
to  obtain  the  finite  difference  form  of  Equation  (6)  as 
(if,  i  -  if )  j(2uX| .  ,  +  uXl)  -  (uXj.,  -  uxj) 


“  (ri>i ,  i  -  ipi) 


(2uPl , ,  +  uPl)  -  (u,,..,  -  u,„) 


r, .  i  +  » , 
ri>, 

rh.i  +  rh 


(H) 

At  the  body  surface  ri  =  r,,,.  This  condition  can  be  used  to  solve 
for  rPl  _  j  step  by  step  from  the  body  surface  outward,  e.g. 


r», 


l‘i .  l 


C  .  V  C-  -  4BD 
2B 


(12) 


u.,  K  f  dtHr)  dr 

—  <=-—/  — —  = - —  •  lit 

V  2  dr  r  tan/?,(r) 

n, 

where  (1  =  r'2nR,,V,  tan/?,  =  (u,.  +  uJ/fQrR,.  -  U|).  K  is^  the 
number  of  blades,  Q  is  the  propulsor  angular _yeloeity,  r  = 
r  K, „  r.  is  the  dimensionless  hub  radius,  and  u.,  and  u*  are  the 
axial  and  tangential  propulsor  induced  velocities  at  lifting 
line,  respectively. 

(cl  A  new  effective  velocity  profile  can  be  obtained  bv 
using  the  new  lifting-line  solution  of  u;l  in  Equation  (9). 

id)  Steps  (b)  and  (cl  are  repeated  until  the  values  of  G 
and  u  converge  to  within  some  specified  tolerance.  Two 
iterations  have  been  found  to  be  sufficient  for  preliminary 
design. 

(e)  If  a  final  lifting  surface  design2"  -2  -11  is  required, 
the  final  solution  of  circumferential-average  propulsor-induced 
axial  velocity,  u„,  at  the  propeller  plane  can  be  calculated  by 
using  a  field  point  velocity  program2"  which  includes  lifting 
surface  corrections  and  thickness  distributions.  The  final  ef¬ 
fective  velocity  profile  is  obtained  by  using  the  final  lifting- 
surface  solution  of  u„  in  Equation  (9). 

For  calculation  of  the  performance  of  a  propulsor  of 
given  geometry  and  rpm.  a  computational  method-’4  is 
available  to  predict  (I  and  tan/?,  for  estimated  values  of  u,.  A 
new  effective  velocity  profile  can  lie  calculated  using  the  new 
induced  velocity  u„  computed  from  the  field-point  velocity 
program.1" 

2  3  Total  Velocities  I  pstream  of  Propulsor 


The  total  velocities  at  the  propulsor  plane  are  rather 
difficult  to  measure.  Therefore,  in  order  to  evaluate  the  pres¬ 
ent  theory  a  comparison  of  total  velocity  u,,  between  theory 
and  exjieriment  must  be  made  immediately  upstream  of  the 
propulsor.  The  total  velocity  is  u,,  =  u,  +  u,,  where  u*  is  solv¬ 
ed  by  Equation  (9)  and  the  value  of  the  propulsor-induced 
velocity  u.,  must  tie  calculated  at  the  position  where  the  total 
velocities  are  measured.  The  new  radial  position  for  Up  can  be 
obtained  from  Equation  (6).  The  nominal  velocity  ux  can  be 
approximated  very  accurately  over  a  small  increment  of 
radius  dr  as  a  linear  function  of  r.  Although  the  velocity  at 
the  wall  is  zero,  the  velocity  profile  in  the  present  approxima¬ 


where  B  =  2u,,, ,  ,  +  uP| 

r  =  -  r!',<ui', .  i  -  %> 

1)  =  -  rf,,(uPl ,  |  -  2uPl)  -  F 

F  =  (r2.  ,  -  rf)(2ux, . ,  +  uX|)  -  r,(ux, .  ,  -  uXl) 

(r, .  i  r.) 

Both  the  effective  velocity  u,.,  and  the  new  radial  position 
rPl  .  |  for  the  total  velocity  u,., ,  ,  =  u.., ,  ,  +  u„, ,  ,  are  the  simple 
solutions  of  the  quadratic  Equations  (9  and  12).  The  effective 
velocity  is  solved  (independent  of  rp)  inward  from  a  large  radial 
position  where  the  propulsor-induced  velocity  is  zero.  Once  the 
values  of  u,  are  obtained,  the  new  radial  position  rIh , ,  can  be 
solved  from  the  body  surface  outward.  These  two  solutions  are 
decoupled  without  any  need  for  employing  an  iterative  numeri¬ 
cal  procedure.  These  solutions  can  readily  be  incorporated  in 
most  propulsor  design  practices. 

3.  EXPERIMENT 

3.1  Two  Propulsors  in  Five  Nominal  Wakes 

The  experimental  data  used  to  evaluate  the  present 
theory  are  derived  from  two  propulsors  operating  in  five  dif¬ 
ferent  nominal  wakes.  The  geometrical  characteristics  of  the 
propulsor  used  at  DTNSRDC,  denoted  Propulsor  T,  are  given 
in  Reference  1  and  those  of  the  propulsor  used  at  Nagasaki 
Experimental  Tank  by  Nagamatsu  and  Tokunaga.  denoted 
Propulsor  J,  are  given  in  Reference  2. 

Three  axisymmetric  Ixidies,  designated  as  Models  1,  C, 
and  D,  were  used  to  generate  thick  stern  boundary  layers  for 
Propulsor  T.  The  offsets  of  the  three  models  are  given  in 
Tables  1  through  3.  The  propulsor  was  located  at  x/L  = 

0.983,  where  x  is  the  axial  distance  from  the  nose  and  L  is 
the  total  body  length.  The  ratio  of  the  propulsor  radius  to  the 
maximum  hull  radius,  Rp/r,,,**.  is  0.545  for  Model  1  and  0.484 
for  Models  C  and  D.  The  maximum  radius  of  the  model  (rnu4X) 
is  13.97  cm  for  Model  1  and  31.18  cm  for  Mixiels  C  and  I). 

The  three  models  were  constructed  of  molded  fiberglass;  spe 
cified  profile  tolerances  were  held  to  less  than  t  0.4  mm.  all 
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imperfections  were  removed,  meridians  were  faired,  and  tne 
fiberglass  was  polished  to  a  0.64  mi  -run  nus  surface  finish. 

One  axisymmetrd  body  was  constructed  by  Nagamatsu 
an1  T  duiiaga-  for  their  investigation.  The  total  length  of  the 
Dody  was  3.085  m  and  the  maximum  radius  was  10.66  cm. 

The  propulsor  was  located  at  2.57  propeller  diameters  down¬ 
stream  of  the  after  perpendicular.  A  wire  mesh  was  wrapped 
around  the  aft  portion  of  the  parallel  middle  body.  The  nom¬ 
inal  wakes  without  and  with  the  wire  mesh  will  be  designated 
as  wakes  A  and  B,  respectively.  The  value  of  Rp/rmiuc  was 
0.610. 

The  measured  distributions  of  the  nominal  axial  veloc¬ 
ities  at  the  propulsor  plane  for  the  five  wakes  are  shown  in 
Figure  2.  There  is  considerable  variation  in  the  nominal  veloc¬ 
ity  distributions  among  the  five  'cakes  selected  for  the  pres¬ 
ent  investigation. 
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Figure  2  -  Measured  Nominal  Axial  Velocity  Profiles 
of  the  Five  Nominal  Wakes 

3  2  Exjienmental  Technique 

The  experimental  investigation  of  Models  1,  t\  and  I) 
was  conducted  in  the  DTNSHIX'  anecboic  wind  tunnel.  The 
wind  tunnel  has  a  2  438-  by  2.438-m  test  section  with  a  maxi¬ 
mum  air  speed  of  61  m/s.  The  model  was  supported  by  two 
streamlined  struts  separated  by  roughly  one-third  of  the 
model  length  The  upstream  strut  had  a  15-cm  chord  and  the 
downstream  strut  had  a  3-cm  chord.  The  distrubances  gener- 
ated  by  the  supporting  struts  were  within  the  region  below  the 
horizontal  centerplane.  All  the  measurements  of  the  nominal 
and  the  total  velocities  were  made  in  the  vertical  centerplane 
along  the  upper  meridian  and  there  was  little  extraneous  ef¬ 
fect  from  the  supporting  struts.  Each  stern  of  the  model  pro¬ 
truded  from  the  closed-jet  working  section  of  the  wind  tunnel 


into  the  anechoic  chamber  (6.4  x  6.4  x  6.4  m)  located  up¬ 
stream  of  the  diffuser.  Propulsor  T  was  driven  by  a  9-kW 
high-speed  motor  mounted  inside  the  stem  of  the  model.  The 
rotational  speed  of  the  propulsor  shaft  was  measured  by  a 
magnetic  pick  o. 

Time-average  nominal  velocities  at  the  propulsor  plane 
and  at  a  plane  immediately  upstream  of  the  propulsor  plane 
in  the  absence  of  the  propulsor  were  measured  by  a  Laser 
Doppler  Velocimeter  (LDV)  and  checked  by  a  TSI,  Inc.  Model 
1241  “X”  wire  monitored  by  a  two-channel  TSI  Model  1050-1 
hot-wire  anemometer."’  The  axial  velocity  profiles  inside  the 
stern  boundary  layer  measured  by  the  LDV  and  by  the  hm 
wire  anemometer  were  found  to  agree  within  two  percent  of 
the  free-stream  velocity.  The  total  velocities  immediately  up¬ 
stream  of  the  propulsor  were  measured  by  the  LDV  with  tin 
propulsor  in  operation.  Some  of  the  data  have  already  in-on 
reported  in  Reference  1.  Most  data  presented  in  this  investi¬ 
gation  were  measured  after  the  publication  of  Reference  1 
although  the  same  Ll)V  measurement  techniques  were  used. 
The  overall  accuracy  of  the  measured  axial  velocities  by  tin 
LDV  is  almut  two  percent  of  the  free-stream  velocity. 

In  the  investigation  of  Nagamatsu  and  Tokimaga.-  the 
model  was  also  supported  by  two  struts  and  propulsor  .1  un¬ 
mounted  on  and  driven  by  the  dynamometer  used  for  open 
water  propulsion  testing  The  experiments  were  conducted  m 
a  towing  tank.  The  propulsor  plane  was  located  at  5.13  K; 
downstream  of  the  after  perpendicular  of  the  model;  the 
model  centerline  depth  of  submergence  was  4.6  R, ,  where  It., 
is  the  radius  of  the  propulsor.  A  five-hole  pitot  probe  of  Nl'l. 
type  was  used  to  measure  the  axial  and  the  radial  velocities 
with  and  without  the  propulsor  in  operation.  The  overall  ac 
curacy  of  this  type  of  pitot  probe  for  measuring  velocity  is 
about  three  percent  of  the  free-stream  velocity.  Standard  self 
propulsion  and  open  water  experiments  were  conducted  to 
determine  the  Taylor  wake  fraction. 

Standard  self-propulsion  experiments  were  also  con¬ 
ducted  in  the  DTNSRDC  model  basin  using  a  large  model;  ihe 
centerline  depth  of  submergence  was  1 7  Ft,  to  avoid  the  cf 
feet  of  the  free  surface.  The  towing  strut  had  a  2<i-cfn  chord 
and  was  located  at  40-percent  of  body  length  downstream 
front  the  nose.  The  conditions  of  the  self-propulsion  expert 
ments  for  the  two  propulsors  operating  in  the  five  nominal 
wakes  are  tabulated  in  Table  4.  to  lie  discussed  below 

4.  RESFLTS  AND  DISCISSIONS 

4.1  Effective  Wake 

The  two  propulsors  used  in  the  self  propulsion  experi¬ 
ments  (Table  4)  were  stock  propulsors  and  were  not  operated 
at  their  design  conditions.  A  propulsor  performance  predic¬ 
tion  computer  program-1  was  used  to  compute  the  values  of 
non-dimensional  circulation  tl  and  hydrodynamic  pitch  angle 
tan/3,  for  the  estimated  values  of  u,..  The  final  values  of  tan  fi, 
were  scaled  up  or  down  by  the  ratio  of  the  measured  value  of 
Kt  to  the  computed  value  of  Kt-  The  values  of  (.1  and  modi¬ 
fied  tan/?,  were  then  used  to  compute  the  propulsor-indueed 
velocities  by  a  field-point  velocity  computer  program-"  with  a 
lifting-surface  option.  The  new  values  of  u„  were  then  used  to 
compute  a  second  estimate  of  the  effective  velocity  u,.  b,  us¬ 
ing  Equation  (9).  Three  iterations  were  sufficient  for  the  com¬ 
puted  values  of  u,.  to  converge  to  within  0.2  percent  accuracy. 
Table  4  shows  the  computed  values  of  Ky.  1i>.  the  volume 
mean  effective  and  nominal  velocities.  (udc  and  tu,  V-  ami  ’be 
effective  velocity  at  the  0.7  propulsor  radius.  The  notation 
used  in  Table  4  is  as  follows: 
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The  exjierimontally  determined  value  of  the  Taylor  wake  frac¬ 
tion  W|  is  also  shown  in  Table  4.  It  is  im|>ortant  to  note  that 
the  measured  values  (l-w-j)  agree  to  within  ±  0.03  with  the 
computed  values  of  the  volume-mean  effective  velocity  ratio, 
|u,.|^/V.  and  effective  velocity  ratio  at  0.7  propulsor  radius, 
|u,p  ;/V,  for  all  twelve  propulsion  conditions.  The  overall  ac¬ 
curacy  of  the  propulsion  experiments  is  about  ±  2%.  There  is 
no  reason  that  the  value  of  (1-w-r)  should  be  exactly  equal  to 
(u.  l^/V  However,  the  general  agreement  between  the 
measured  values  of  (l-wT)  and  the  computed  values  of  |u,.|yV 
lor  a  wide  variety  of  nominal  wakes  and  propulsion  loading 
conditions  does  provide  an  indirect  confirmation  of  the  pres¬ 
ent  theory. 

Typical  values  of  the  computed  effective  and  total 
velocities  at  the  plane  of  the  propulsor  are  shown  in  Figure  3 
and  listed  in  Table  f>.  The  final  computed  values  of  (I,  tan/3,, 
and  u,.  are  listed  in  Table  II  for  the  conditions  used  in  this  in 
vestigation.  This  information  is  sufficient  for  other  in¬ 
vestigators  to  perform  independent  calculations. 


1.2  Total  Velocities  Upstream  of  I’ropulsor 


With  the  final  computed  values  of  tl,  tan/3,  and  u,.,  the 
field-point  velocity  computer  program-’"  with  the  lifting- 
surface  option  was  used  to  compute  the  induced  velocity  u„ 
upstream  of  the  propulsor.  The  total  velocity  n,,  and  the 
radial  position  r,,  were  then  computed  using  Kquations  (9)  and 
(1 11.  The  detailed  results  of  the  computation  and  the  meas¬ 
ured  data  are  tabulated  in  Table  7.  The  comparisons  iielween 
the  measured  total  velocity  profiles  and  the  computed  total 
velocity  profiles  are  shown  in  Figures  4  through  9.  The  agree¬ 
ment  between  the  measured  and  computed  values  of  u,,  is  ex 
eellent  for  all  cases.  Thus,  it  may  lie  concluded  that  the  pres 
ent  inviseid  approximation  provides  a  good  representation  of 
the  complex  hydrodynamic  interaction  Iielween  the  propulsor 
and  the  wake  and  can  be  used  with  confidence  to  calculate 
the  effective  velocity  profile  for  a  propulsor  in  an  axisym- 
metric  wake.  The  ratio  of  the  measured  total  and  nominal 
velocities  shown  in  Figures  1  through  9  is  found  to  be  largest 
near  the  propulsor  hub  and  to  decrease  toward  the  propulsor 
tip.  Hence,  the  computed  ratio  of  effective  and  nominal 
velocities  is  largest  at  the  propulsor  hub  and  decreases  to 
about  l.i>  at  the  propulsor  tip,  which  is  quite  different  from 
methods  which  assume  that  the  effective  wake  profile  is  a 
constant  multiple  of  the  nominal  wake  profile. 

f».  CONCLUSION 

In  this  paper,  we  have  summarized  recent  experi¬ 
mental  and  theoretical  investigations  of  the  effective  wake  for 
the  axisymmetric  propulsor/stern  boundary  layer  interaction 
problem,  A  comprehensive  set  id’  experimental  data  in  tabu¬ 
lated  form  is  presented  lor  two  propulsors  o|H>rating  in  five 
quite  different  nominal  wakes.  The  data  are  compared  with 
calculations  based  on  an  inviseid  propulsor'stern  boundary 
layer  interaction  theory  and  provide  valuable  insights  into  the 
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manner  in  which  the  propulsor  interacts  with  a  thick  stern 
boundary  layer. 

The  numerical  procedures  for  computing  the  effective 
wake  and  the  total  velocity  upstream  of  the  propulsor  have 
been  simplified  to  the  solution  of  two  quadratic  algebraic 
equations.  The  procedure  can  readily  be  incorporated  into  any 
propulsor  design  and  performance  prediction  practice. 

It  is  shown  that  for  twelve  self-propulsion  conditions 
covering  a  wide  range  of  propulsion  thrust  loading  coeffici¬ 
ents,  the  computed  volume-mean  effective  wakes  agree  well 
with  the  measured  Taylor  wake  fractions  based  on  thrust 
identity.  The  present  inviseid  propulsor/stern  boundary  layer 
interaction  theory  predicts  very  well  the  measured  total 
velocity  profile  immediately  upstream  of  the  operating  pro¬ 
pulsor  The  computed  ratio  of  the  effective  and  nominal  veloc¬ 
ities  is  largest  near  the  propulsor  hub  and  decreases  to  alviut 
1.0  at  the  propulsor  tip.  This  is  contrary  to  some  older  meth¬ 
ods  which  assume  that  the  effective  wake  profile  is  a  constant 
multiple  of  the  nominal  wake  profile. 
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Figure  7  -  Measured  and  Computed  Axial  Velocity  f'rofiles 
Immediately  Upstream  of  Propulsor  (x/I)(,  -  -0.227) 
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Figure  8  —  Measured  and  Computed  Axial  Velocity  Profiles 
Immediately  Upstream  of  Propulsor  (x/Dp  -  -0.239) 
With  and  Without  Propulsor  T  Operated  in  Wake  C 
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FRACTION  AT  VARIOUS  PROPULSION  CONDITIONS  IN  THE  FIVE  WAKES 
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TABU  •  -  TOTAL  ANO  EFFECTIVE  AXIAL  VELOCITY  PROFILES  COMPUTED  FROM  NOMINAL  AND 
PROPULSOR  INDUCED  AXIAL  VELOCITIES  AT  THE  LOCATION  OF  PROPULSOR 


TABLE  1  -  COMPUTED  RESULTS  OF  PROPULSOR/WAKE  INTERACTION 
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TABLE  7  -  MEASURED  AND  COMPUTED  VELOCITY  PROFILES  IMMEDIATELY  UPSTREAM  OF  PROPULSOR 
WITH  AND  WITHOUT  AN  OPERATING  PROPULSOR 
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